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Abstract 

Oh . 

^ We investigate the dynamics of the Fisher equation for the spreading of micro-organisms in one dimenison subject to both turbulent 
• i-H convection and diffusion. We show that for strong enough turbulence, bacteria , for example, track in a quasilocalized fashion 
i-^ (with remakably long persistance times) sinks in the turbulent field. An important consequence is a large reduction in the carrying 

^Htapacity of the fluid medium. We determine analytically the regimes where this quasi-localized behavior occurs and test our 
I— I predictions by numerical simulations. 
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The spreading of bacterial colonies at very low Reynolds 
numbers on a Petri dish can often be described [Tj by the 
Fisher equation [2], i.e. 

\ dtc = Ddl^c + fic — bc^ , (1) 

where c(x, t) is a continuous variable describing the con- 
centration of micro-organisms, D is the diffusion coefficient 
and /i the growth rate. 

' In the last few years, a number of theoretical and ex- 
perimental studies [3] , [3| , [S] , [5] , [7] have been performed 
to understand the spreading and extinction of a popula- 
tion in an inhomogeneous environment. In this paper we 
study a particular time-dependent inhomogeneous enviro- 
ment, namely the case of the field c{x,t) subject to both 
convection and diffusion and satisfying the equation: 

dtc + div{Uc) = DV^c + ^ic-bc^ (2) 

where U{x, t) is a turbulent velocity field. Upon specializing 
to one dimension, we have 

dtc + d^iUc) Ddlc + ^ic-bc^ (3) 

Equation ([3]) is relevant for the case of compressible flows, 
where d^U ^ 0, and for the case when the field c{x,t) 
describes the population of inertial particles or biological 
species. For inertial particles, it is known [9] that for large 
Stokes number , i.e. the ratio between the characteristic 
particle response time and the smallest time scale due to 
the hydrodynamic viscosity, the flow advecting c{x, t) is 



effectively compressible, even if the particles move in an 
incompressible fluid. Let us remark that the case of com- 
pressible turbulence is also relevant in many astrophysical 
applications where ([2|) is used as a simpliflcd prototype of 
combustion dynamics. By suitable rescaling of c{x.t), we 
can always set 5 = 1. In the following, unless stated oth- 
erwise, we shall assume 5=1 whenever n ^ and 6 = 
for ^ = 0. For a treatment of equation ([3]) with a spatially 
uniform but time-dependent random velocity, see [8] 

The Fisher equation has travelling front solutions that 
propagate with velocity vp ^ (^m)^^^ E], [ini- In Fig. ^ 
we show a numerical solution of Eq. ([T]) with D ~ 0.005, 
^ = 1 obtained by numerical integration on a space do- 
main of size L = 1 with periodic boundary conditions. The 
figure shows the space-time behaviour of c(x,t), the color 
code representing the curves c(x, t) = const. With initial 
condition c{x, t = 0) nonzero on only a few grid points cen- 
tered at X = L/2, c{x, t) spreads with a velocity vp ^ 0.07 
and, after a time L/vp ~ 4 reaches the boundary. 

A striking result, which motivated our investigation, is 
displayed in Fig. ([2]), showing the numerical solutions of Eq. 
([3]) for a relatively "strong" turbulent flow, where the av- 
erage convection velocity vanishes and " strong turbulence" 
means high Reynolds number ( a more precise definition of 
the Reynolds number and specification of the velocity field 
is given in the following sections). From the figure we see no 
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Fig. 1. Numerical simulation of cq. [T]with fi = 1, D = 0.005 and with 
periodic boundary conditions. The initial conditions are c{x,t) = 
everywhere expect for few grid points near L/2 = 0.5 where c = 1. 
The horizontal axis represents time while the vertical axis is space. 
The colors display different contour levels of c(x,t). 

trace of a propagating front: instead, a well-localized pat- 
tern of c(x, t) forms and stays more or less in a stationary 
position. 

For us, Fig. ([2|) shows a counter intuitive result. One naive 
expectation might be that turbulence enhances mixing. The 
mixing effect due to turbulence is usually parametrized in 
the literature [TT] by assuming an effective (eddy) diffusion 
coefficient D^ff 3> D. As a consequence, one naive guess 
for Eq. ([3]) is that the spreading of an initial population 
is qualitatively similar to the travelling Fisher wave with 
a more diffuse interface of width ^Jj^fjlJi. As wc have 
seen, this naive prediction is wrong for strong enough tur- 
bulence: the solution of equation shows remarkable lo- 
calized features which are preserved on time scales longer 
than the characteristic growth time or even the Fisher 
wave propagation time L/vp. An important consequence of 
the localization effect is that the global " mass" (of growing 
microorganisms, say) , Z = J dxc{x, t), behaves differently 
with and without turbulence. In Fig. ([3]), we show Z{t): the 
curve with red circles refers to the conditions shown in Fig. 
Il])), while the curve with green triangles to Fig. 

The behavior of Z for the Fisher equation without tur- 
bulence is a familiar S'-shaped curve that reaches the max- 
imum Z = 1 on a time scale L/vp. On the other hand, the 
effect of turbulence (because of localization) on the Fishe 
equation dynamics reduces significantly Z almost by one 
order of magnitude. 

With biological applications in mind, it is important 
to determine conditions such that the spatial distribution 
of microbial organisms and the carrying capacity of the 
medium are significantly altered by convective turbulence. 
Within the framework of the Fisher equation, localization 
effect has been studied for a constant convection velocity 
and quenched time-independent spatial dependence in the 
growth rate /i [6], [7], [H], [13] • In our case, localization, 
when it happens, is a time-dependent feature and depends 
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Fig. 2. Same parameters and initial condition as in Fig. |(TJ for 
equation l[3} with a "strong turbulent" flow u advecting c{x,t). 
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Fig. 3. The behavior in time of the total "mass" Z(t) = J dxc{x, t). 
The red circles show the function Z for the case of Fig. JTJ, i.e. a 
Fisher wave with no turbulence. The green triangles show Z for the 
case of Fig. ((2]l when a strong turbulent flows is advecting c{x,t). 

on the statistical properties of the compressible turbulent 
flows. As discussed in detail below, a better term for the 
phenomenon we study here might be " quasilocalization" , in 
the sense that (I) spatial localization of the growing popula- 
tion sometimes occurs at more than one location; (2) these 
spatial locations drift slowly about and (3) localization is 
intermittent in time, as localized populations collapse and 
then reform elsewhere. For these reasons, the quasilocaliza- 
tion studied here is not quite the same phenomena as the 
Anderson localization of electrons in a disordered poten- 
tial studied in [TB]. Nevertheless, the similarities are suffi- 
ciently strong that we shall use the terms "quasilocaliza- 
tion" and "localization" interchangeably in this paper. It 
is worth noting that the localized "boom and bust" pop- 
ulation cycles studied here may significantly effect "gene 
surfing" [TJ] at the edge of a growing population, i.e. by 
changing the probability of gene mutation and fixation in 
the population. 



For the case of bacterial populations subject to both tur- 
bulence and convections due to, say, an external force such 
as sedimentation under the action of gravity, we may think 
that the turbulent velocity can be decomposed into a con- 
stant "wind" Mo and a turbulent fluctuation u{x,t) with 
zero mean value, U{x,t) = uq + u{x,t). We find that the 
localization shown in Fig. ^ can be significantly changed 
for large enough background convection uq. 

We would like to understand why and how uq 7^ can 
change the statistical properties of c{x, t) in the presence 
of a random convecting velocity field. We wish to under- 
stand, in particular, whether c{x, t) spreads or localizes as 
a function of parameters such as the turbulence intensity 
and the mean "wind" speed uq. 

Our results are based on a number of numerical simula- 
tions of Eq. Q performed using a particular model for the 
fluctuating velocity field t). In Sec. 1 we introduce the 
model and we describe some details of the numerical simu- 
lations. In Sec. 2 we develop a simple " phenomenological" 
theory of the physics of Eq. ^ based on our present under- 
standing of turbulent dynamics. In Sec. 3 we analyze the 
numerical results when the sedimentation velocity uq ~ Q 
while in Sec. 4 we describe our findings for uq > . Con- 
clusions follow in Sec. 5. 

1. The model 

To completely specify equation ([3]) we must define the 
dynamics of the "turbulent" velocity field U{x, t). For now, 
we set U{x,t) = uq + u{x,t), neglect the uniform part 
uq = and focus on u{x,t). Although we consider a one 
dimensional case, we want to study the statistical proper- 
ties of c{x, t) subjected to turbulent fiuctuations which are 
close to thos generated by the three dimensional Navier- 
Stokes equations. Hence, the statistical properties oiu[x, t) 
should be described characterized by intermittency both 
in space and in time. We build the turbulent field u[x, t) 
by appealing to a simplified shell model of fiuid turbulence 
[T5] . The wavenumber space is divided into shells of scale 
fc„ = 2"~^A:o, n = 1,2,.... For each shell with character- 
istic wavenumber A:„, we describe turbulence by using the 
complex Fourier- like variable M„(t), satisfing the following 
equation of motion: 

+ (1 - (5)fc„_ilt„„iU„_2) + fn ■ (4) 

The model contains one free parameter, (5, and it conserves 
two quadratic invariants (when the force and the dissipa- 
tion terms are absent) for all values of 5. The first is the 
total energy |u„p and the second is X]n(~l)"''''n I"" 
where a = log2(l — 5). In this note we fix 5 = —0.4. For 
this value of 8 the model reproduces intermittency features 
of the real three dimensional Navier Stokes equation with 
surprising good accuracy |15| . Using u„, we can build the 
real one dimensional velocity field u{x,t) as follows: 



u{x,t) = F^Ke'^-"^^ +<e-*'="1, (5) 

n 

where is a free parameter to tune the strength of velocity 
fluctuations (given by u„) relative to other parameters in 
the model (see next section). In all numerical simulations 
we use a forcing function /„ = (e{l + i)/ul)Sn.i, i-G- energy 
is supplied only to the largest scale corresponding to n = 

1. With this choice, the input power in the shell model is 
simply given by 1/2 J2nlKfn+Unfn] = e > i-S- it is constant 
in time. To solve Eqs. ([3]) and ^ we use a finite difference 
scheme with periodic boundary conditions. 

Theses model equations can be studied in detail without 
major computational efforts. One main point of this note 
is to explore the qualitative and quantitative dynamics of 
Eqs. ©,(111) and compare it against the phenomenological 
theory developed in the next section. 

The free parameters of the model are the diffusion con- 
stant D, the size of the periodic Id spatial domain L, the 
growth rate /i, the viscosity f (which fixes the Reynolds 
number Re), the mean constant vclociy ug, the "strength' 
of the turbulence F and finally the power input in the shell 
model, namely e. Note that according to the Kolmogorov 
theory [iT], e ^ u^^^/L where u^ms is the mean square 
velocity. Since Urms ~ F, we obtain that F and e are re- 
lated as e ^ F^. By rescaling of space, we can always put 
L = 1. We fix e = 0.04 and v = 10~^, corresponding to an 
equivalent Re = UrmsL/v ~ 3 x 10^. As we shall see in the 
following, most of our numerical results are independent of 
Re when Re is large enough. In the limit Re — > 00, the sta- 
tistical properties of eq. ([3]) depend on the remaining free 
parameters, D, uq, fi and F. The important combinations 
of these parameters are discussed in the next section. 

2. Theoretical considerations 

We start our analysis by rewriting ([3]) in the form: 

dtc + (uo + w)d^ (c) = Dd^c +ip + g)c - hc^ (6) 

where w = u{x,t) and g{x,t) = -~dxu{x,t). Previous the- 
oretical investigations [5] have shown that for uq = w = 
0, c(x,t) becomes localized in space for time-independent 
"random" forcing g = g{x) (Anderson localization [l6]). 
For Mo large enough, a transition from localized to extended 
solutions has been predicted and observed in previous nu- 
merical and theoretical works [12]. Here, we wish to un- 
derstand whether something resembling localized solutions 
survives in equation ([6]) when both w and g depend on time 
as well as space. 

To motivate our subsequent analysis, consider first the 
case /i = lib = 0. In this limit, Eq. ([3]) is just the Fokker- 
Planck equation describing the probability distribution 
P{x, t) = c{x, t) to find a particle in the range (x, x + dx) 
at time t, whose dynamics is given by the stochastic differ- 
ential equation: 

dx / 

— ^u{x,t) + V2Di^{t) (7) 
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where r]{t) is a white noise with {ri{t)ij{t')) = 6{t — t'). 
Let us assume for the moment that u(x,t) — u{x) is time 
independent. Then, the stationary solution of ([3]) is given 

by 

P{x,t) = A-^exp[-<S>{x)/D] (8) 

where A is a normahzation constant and = —u{x). 
It follows that P{x,t) = P{x) is strongly peaked near the 
points Xi where $ has a local minimum, i.e. u{xi) = and 
—d^x^ = dxu{x)\x=xi < 0. Let us now consider the be- 
haviour of P{x) near one particular point xq where u{xo) = 
0. For X close to xq wc can write: 

-To{x ~ xo) + V2DTj{t) (9) 



dt 

where Fq = —dxu{x)\x=xo- Equation ^ is the Langevin 
equation for an overdamped harmonic oscillator, and tells 
us that P is spread around xq with a characteristic " local- 
izaiont length" of order ^/ = ^JdJTq. On the other hand, 
wc can identify Fq with F, a typical gradient of the turbu- 
lent velocity field u. In a turbulent flow, the velocity field 
is correlated over spatial scale of order i'*/F where v1/2 is 
the average kinetic energy of the flow. For P to be local- 
ized near xq, despite spatial variation in the turbulent field, 
we must require that the localization length ^; should be 
smaller than the turbulent correlation scale f*/F, i.e. 



DT 



> 2 



(10) 



Condition (|10p can be easily understood by considering 
the simple case of a periodic velocity field u, i.e. u = 
u*cos(xi;*/F). In this case, condition (fTO)) states that D 
should be small enough for the probability P not to spread 
over all the minima of u. For small D or cquivalently for 
large wJ/F, the solution will be localized near the minima 
of u, at least for the case of a frozen turbulent velocity field 
u{x). 

The above analysis can be extended for velocity field 
u(x, t) that depend on both space and time. The crucial 
observation is that, close to the minima Xi of <f>(a;,t) = 
— J dxu{x, t), we should have u{xi, t) ~ 0. Thus, although 
u is a time dependent function, sharp peaks in P{x, t) move 
quite slowly, simply because u{x, i) ^ near the maximum 
of P{x,t). One can consider a Lagrangian path x{t) such 
that x{0) ~ zq, where xq is one particular point where 
u(a;o,0) = and ^^^^(a;, O)!^;^^:^ < 0. From direct numer- 
ical simulation of Lagrangian particles in fully developed 
turbulence, we know that the acceleration of Lagrangian 
particles is a strongly intermittent quantitiy, i.e. it is small 
most of the time with large (intermittent) bursts. Thus, 
we expect that the localized solution of P follows x{t) for 
quite long times except for intermittent bursts in the tur- 
bulent flow. During such bursts, the position where u ~ 
changes abruptly, i.e. almost discontinuosly from one point, 
say x{t), to another point x{t + St). During the short time 
interval (5i, P will drift and spread, eventually reforming to 



become localized again near x{t + 5t). The above discus- 
sion suggests that the probability P{x, t) will be localized 
most of the time in the Lagrangian frame, except for short 
time intervals 5t during an intermittent burst. 



We now revisit the condition (|T0|) . For the case of a time- 
dependent velocity held u, we estimate F as the character- 
istic gradient of the velocity field, i.e. 

T^{{dxufY'^ 
where (..) stands for a time average. Now, v1 should be 
considered as the mean kinetic energy of the turbulent fluc- 
tuations. In our model, both and F are proportional to 
the strength of the velocity fluctuations. Thus, we can 
rewrite the localization criteria (1101) in the form: 
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> 2 



(11) 



where and F are computed for F = 1. We conclude that 
for small values of F, P{x,t) is spread out, while for large 
F, P should be a localized or sharply peaked function of x 
most of the time. An abrupt transition, or at least a sharp 
crossover, from extended to sharply peaked functions P, 
should be observed for increasingf F. 

It is relatively simple to extend the above analysis for a 
non zero growth rate fi > 0. The requirement PH)) is now 
only a necessary condition to observe localization in c. For 
/.J > we must also require that the characteristic gradient 
on scale must be larger than i.e. the effect of turbulence 
should act on a time scale smaller than 1/fi. We estimate 
the gradient on scale as Sv{^i)/^i, where Sv{^i) is the 
characteristic velocity difference on scale ^; . We invoke the 
Kolmogorov theory, and set Sv{S,i) = v^{(^i/ L)^/^ to obtain: 

M<^-^-W^)^/^ (12) 

In (fT^ , we interpret F as the characteristic velocity gradi- 
ent of the turbulent flow. Because v^, ^ F and F ~ i^, it 
follows that the r.h.s of (fT2|) goes as F^^^. Note also that 
^"(6)/6 < F on the average, which leads to the inequality: 

fi<T (13) 

From (IlOl and (Ull) we also find 



Dfi 



> 2 



(14) 



a second necessary condition. Once again, we sec that lo- 
calization in a Lagrangian frame should be expected for 
strong enough turbulence. 

One may wonder whether a non zero growth rate fi can 
change our previous conclusions about the temporal be- 
havior, and in particular about its effect on the dynamics 
of the Lagrangian points where u{x, t) = 0. Consider the 
solution of ([3]) at time t, allow for a spatial domain of size 
L, and introduce the average position 
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where Z{t) = dxc{x,t). Upon assuming for simplicity 
a single localized solution, we can think of Xm just as the 
position where most of the bacterial concentration c{x, t) is 
localized. Using Eq. ^ , we can compute the time derivative 
Vm(t) = dxm/dt. After a short computation, we obtain: 

(xm. — x)P{x,t)'^ + / u{x,t)P{x,t)dx 
Jq 

(16) 

where P{x,t) = c{x,t)/Z{t) and Z{t) = c{x,t)dx. Note 
that Vm is independent of ^. Moreover, when c is localized 
near Xm, both terms on the r.h.s. of (fT6|) are close to zero. 
Thus, Vm can be significantly different from zero only if 
c is no longer localized and the first integral on the r.h.s 
becomes relevant. We can now understand the effect of the 
non linear term in when c{x,t) is localized, the non 
linear term is almost irrelevant simply because is close 
to 0. On the other hand, when c{x,t) is extended the non 
linear term drives the system to the state c = 1 which is 
an exact solution in the absence of turbulent convection 
u{x,t) ~ 0. 

We now allow a non zero mean flow uq ^ 0. As before, 
we first set ^ = 6 = and consider a time -independent 
velocity field u{x). Since the solution c{x, t) of ^ can still 
be interpreted as be the probability to find a particle in the 
interval [x, x + dx\ at time t, we can rewrite ([9]) for the case 
uo > as follows: 

— = -ro(a; - a-o) + + V2l5T^{t) (17) 

The solution of ([3]) is localized near the point xi — xq + 
uq/Tq. Thus for small uq or large Fq there is no major 
change in the arguments leading to (fTO|) . In general, we 
expect that P{x) will be localized near x ~ I'o, provided 
the length = \xi — xo\ = uo/Fo is smaller than ^j, i.e. 
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When p8|) is satisfied, then our previous analysis on local- 
ized solutions for both /.j = and /i 7^ is still valid. Let 
us note that by combing and p8|) we obtain 



"0 



< 2 



(19) 



as a condition for localization, obtained in the study of lo- 
calized/extended transition for steady flowsin an Eulcrian 
context [I], [12]. Here we remark that in a turbulent flow, 
Eq. ([T9|) is only a necessary condition, because (fTO|l . (fT2|) and 
must all also be satisfied for c to show quasi-localized 
states. 

To study the change in the spatial behaviour of P as a 
function of time, we need a measure of the degree of local- 
ization, to look for an some kind of order parameter. Al- 
though there may be a number of valuable solutions, an ef- 
ficient measure should be related to the " order" /" disorder" 
features of c{x, t), where "order" means quasi-localized and 
"disorder" extended. As pointed out in the introduction. 
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Fig. 4. The behavior in time of S{i) . The red circles show the 
function S{t) for the case of Fig. JTJl, i.e. a Fisher wave with no 
turbulence. The blue triangles show Sit) for the case of Fig. ((2j 
when a "strong" turbulent flows is advecting c(x,t). 

the total "mass" of the organisms Z{t) = dxc{x,t) is 
strongly affected by a strongly peaked (or quasi-localized) 
c(x, t), as opposed to a more extended concentration field. 
However, a more illuminating quantity, easily studied in 
simulations, is 



S{t) = - /" dxP{ 
Jo 



x,t)log{P{x,t)) 



(20) 



where P{x,t) = c{x,t)/Z. Localized solutions of Eq. ([3]) 
correspond to small values of this entropy-like quantity 
while extended solutions correspond to large values of S', 
which can be interpreted as the information contained in 
the probability distribution P{x,t) at time t. In our nu- 
merical simulations, we consider a discretized form of ([20|). 
namely: 



Sit) = - E 

i=l,N 



±^logC±^) (21) 



where Xi are now the N grid points used to discretize ([3]) , 
c{xi^t) is the solution of ([3|) in Xi at time t and Z[t) = 
E» c(x,;,t). 

To understand how well S{t) describes whether c{x, t) is 
localized or extended, we consider the cases discussed in 
the Introduction in Fig.s ([T|) and ([2|). The numerical com- 
putations were done with F = for Fig. ([T|) and F ~ 0.9 
for Fig. i.e. no turbulence and "strong" turbulence (the 
attribute "strong" refers to the conditions (fTO|) and (fT2|) . 
In Fig. ^ we show S{t) corresponding to the two simula- 
tions, namely F = (red circles) and F = 0.9 (blue tri- 
angles). The initial condition is the same for both simula- 
tions c(.T,0) = exp[—{x — L/2)^/0.05], i.e. a rather local- 
ized starting point. It is quite clear, from inspecting Fig. 
([ll, that S{t) is a rather good indicator to detect whether 
c(x, t) remains localized or becomes extended. While for 

= (a quiescent fluid), S reaches its maximum value( 
S = 9/05(2.) for 512 grid points) at t = 6. (corresponding 



to uniform concentration c{x,t) ~ 1), for F = 0.9, S is al- 
ways close to its initial value S' ~ 4., indicating that c{x, t) 
is localized, in agreement with Fig. ([2]). 

Let us summarize our findings: when subjected to tur- 
bulence, we expect c{x,t) to be "localized", i.e. strongly 
peaked, most of the time for large enough F and uq = 
/i = 6 = 0. Upon increasing the growth rate /i, the value 
of F where c{x,t) shows Lagrangian localization should 
increase. Finally, for fixed F and /i we should find a lo- 
calized/extended crossover for large enough values of uq. 
Because our theoretical analysis is based on scaling argu- 
ments, we are not able to fix the critical values for which 
localized/extended transition should occour as a function 
of D,F and uq. However, we expect that the conditions 
(fT0 |) .(|12 |) and (fT8|) capture the scaling properties in the pa- 
rameter space of the model. Finally, we have introduced 
an entropy like quantity S{t) useful for analyzing the time 
dependence of of c{x, t) and for distinguishing between lo- 
calized and extended solutions. In the following section, we 
compare our theoretical analysis against numerical simula- 
tions. 



3. Numerical results for Uq — 

We now discuss numerical results obtained by integrat- 
ing equation Eq. As discussed in Sec. 1, all numerical 
simulations have been done using periodic boundary condi- 
tions. Eq. ([3]) has been discretized on a regular grid oi N = 
512 points. Changing the resolution N, shifting N to N ~ 
1024 or = 128, does not change the results discussed in 
the following. We use the same extended initial condition 
c{x, t) — 1 for all numerical simulations with the few excep- 
tions which discussed in the introduction (the Fisher wave) 
and in the conclusions. For all simulations studied here, the 
diffussion constant D has been kept fixed a.t D = 0.005. 

We first discuss the case fi = b = term and beginby un- 
derstanding how well S{t) describes the localized/extended 
feature of the c{x,t). In Fig. ^ we plot S{t) as a func- 
tion of time for a case with F = 0.5. . The behaviour of 
S{t) is quite chaotic, as expected. In Fig. ([6]) we show the 
functions c{x, t) for two particular times, namely t = 35. 
(lower panel) and t = 60 (upper panel). These two partic- 
ular configurations correspond to extended (t = 35) and 
quasi-localized {t = 60) solutions. The corresponding val- 
ues of S are S = 5.5 for t = 35. and S = 3.5 for t = 60. 
It is quite clear that for small S strong localization charac- 
terizes c{x) while for increasing S the behaviour of c{x) is 
more extended. 

In Fig. ^ we also show (red circles) the instantaneous 
behavior of u{x,t) (multiply by a factor 10 to make the 
figure readable). As one can see, the maximum of c{x,t) 
always corresponds to points where u = 0. 
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Fig. 5. The behavior of S{t) for a numerical simulation of Q for 
/i = 6 = UQ = 0, no saturation term c^{x.t) and F = 0.5. With our 
uniform initial condition, 5(0) = 9Zog(2.). As discussed in Sec. 2, S{t) 
is a reasonable indicator for localized/extended spatial behavior of 
c{x, t). Since the flow is turbulent, S{t) behaves chaotically. However, 
it fluctuates at values lower that S{t = 0) and indicates the degree 
of localization. 



t=60. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



space 

Fig. 6. Numerical simulation of Eq. (jsjl for fi = uo = b = 0, term and 
-F = O.5.. The upper panel shows c{x,t) (blue triangles) and u{x,t) 
(red circles) multiply by 10, at t = 60.. The lower panel shows the 
same quantities at t = 35.. The two time frames have been chosen to 
illustrate localized (upper panel) and more extended (lower panel) 
solutions. Note that the localized solution at i = 60. reaches its 
maximum value for u = and d^ulu—o < 0, as predicted by the 
analysis of sec. 2. 

To understand whether the analysis of Sec. 2 captures 
the main features of the dynamics, we plot in Fig. ([7]), 
for t = 35 (lower panel) and t = 60 (upper panel), the 
quantity P{x,t) as computed from Eq. ([H]), i.e. by using 
the instantaneous velocity field u{x,t). Although there is a 
rather poor agreement between c{x, t) and P{x, t) at t = 35, 
at time t = 60 the P{x^ t) is a rather good approximation 
of c{x, t), i.e. when c{x, t) is localized. The spatial behavior 
of c{x, t) is dictated by the point xq where u — and 
the velocity gradient du\x=xo is large and negative. All the 
above results are in qualitative agreement with our analysis. 
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Fig. 7. Numerical simulation of ^ for fi = b = uq = and _F = O.5.. 
The upper panel shows c{x.t) (red circles) at t = 60. and the behavior 
of P{x, t) (blue triangle) computed using ((Sjl using the instantaneous 
velocity field shown in Fig. In the lower panel we show the same 
quantities {c{x, t) and P{x, t)) at time / = 35. when the solution is 
extended. 

Next we test the condition (fTTj) . which states that local- 
ization should become more pronounced for increasing val- 
ues of _F. To test Eq. (jlip we performed a number of numer- 
ical simulations with long enough time integration to reach 
statistical stationarity. In Fig. ([8]) we show {S) as function 
of F, where (...) means a time average. In the insert of the 
same figure, we show to time dependence of S{t) for two 
different values of F, namely F = 0.4 and F = 1.8. The 
behavior of {S) is decreasing as a function of F, in agree- 
ment with ()lip . The temporal behavior of S, for two indi- 
vidual realizations shown in the insert, reveals that, while 
on the average 5* decreases for increasing F, there are quite 
large oscillations in S, i.e. the system shows both localized 
and extended states during its time evolution. However, for 
large F localization is more pronounced and frequent. On 
the other hand, for small values of F, localization is a " rare" 
event. Overall, the qualitative picture emerging from Fig. 
([8]) is in agreement with Eqs. (fTO|) and (fTTj) 

In the previous section we argued that (fTO|l and (fTTj) 
apply also for a time-dependent function velocity u{x,t). 
The basic idea was that c(x, t) is localized near some point 
xq which slowly changes in time, except for intermittent 
bursts. In sufficiently large systems, localization about mul- 
tiple points is possible as well. During the intermittent 
burst, c{x, t) spreads and after the burst c(x, t) becomes lo- 
calized around a new position xq. We have already shown, 
in Figs. ^ and ([7]), that our argument seems to be in 
agreement with the numerical computations using a time- 
dependent velocity field. To better understand this point, 
we measure Vm defined in Eq. (fTB]) . We expect a small Vm 
during localized epochs when S is small. Each time interval 
when c is localized, should end and start with an intermit- 
tent burst where \ vm \ may become large. Figs ^ illustrates 
the above dynamics. The solid red curve is Vm{t) multi- 
plied by a factor 10 while the blue dotted curve shows S{t). 
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Fig. 8. Time averaged entropy {S) as a function of -F. For large 
F the system fluctuates about small values of S,i.e. c{x, t) becomes 
more localized. In the insert, we show the time behaviour of S{t) for 
two particular values of F, namely F = 0.4 (red curve) and F = 1.8 
(blue curve). Numerical simulations performed for fi = b = uq = 

The numerical simulation is for F = 2., i.e. to a case where 
localization is predominant in the system. Fig. ^ clearly 
shows the "intermittent" bursts in the velocity Vm- The 
stagnation point velocity, punctuated by large positive and 
negative excursions, typically wanders near 0. If we assume 
a single sharp maximum in c(x, t), as in the upper panel of 
Fig. the localized profile c{x, t) does not move or moves 
quite slowly. During an intermittent burst, v„i grows signif- 
icantly while c(x, t) spreads over the space. Soon after the 
intermittent burst (see for instance the snapshot at time 
t = 15 in Fig. ([9|)), the velocity v„i becomes small again and 
the corresponding value of S decreases. Fig. ^ provides a 
concise summary of the dynamics: both localized and ex- 
tended configurations of c(x, t) are observed as a function 
of time. During a era of localization, a bacterial concen- 
tration described by c{x,t) is in a kind of "quasi-frozen" 
configuration. 

Fig. ^ tells us that condition (|10p . which was derived 
initially for a frozen turbulent field m, works as well for 
time dependent turbulent fluctuations. As F increases, the 
system undergoes a sharp crossover and the dynamics of 
c{x, t) slows down in localized configurations. Additional 
features of this transition will be discussed later on when 
we focus on a quantity analogous to the specific heat. 

Finally in Fig. (|10p we show the probability distribution 
P{S), obtained by the numerical simulations, for three dif- 
ferent values of F, namely F = 0.2, 0.8 and F = 2.. As one 
can see, the maximum P{S) is shifted toward small values 
of S for increasing F, as we already know from Fig. ([5]). 
Fig. (fTO|) shows that the fluctuations of 5* about the mean 
are approximately independent of F. 

We now turn our attention to the case /i > . We 
have performed numerical simulations for two growth rates, 
namely /i = 1 and /i = 5. We start by analyzing the results 
for ^ = 1. In Fig. pT|) . we show the behaviour of (S) as 
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Fig. 9. Time dependence of Vm (red curve) and S(t) (blue curve) 
for F = 2.0. The value of Vm is multiplied by 10 to make the figure 
readable. The velocity of the accumulation point for the bacterial 
concentration c{x,t), Vm , is computed using l)16|l. We again set 
/I = uo = b = 0. 
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Fig. 10. Probability distribution P{S) of "entropy" S defined by 
II21I I. obtained by the numerical simulations, for three different values 
of F, namely F = 0.2, 0.8 and F = 2.. 

a function of F, while in the insert we show the probabil- 
ity distribution P{S) for three values of S. Upon compar- 
ing with Fig. (fTTj) against Fig.s ^ and (fTU|) . we see that a 
nonzero growth rate /i = 1 does not change the qualitative 
behavior of the system, in agreement with our theoretical 
discussions in the previous section. 

It is interesting to look at the time averaged bacterial 
mass {Z) as a function of F. In Fig. we show {Z) and 
(S) I Smax 8-8 a function of F. For large F, when localization 
dominates the behavior of c{x,t), (Z) is quite small, order 
0.1 of its maximum value, i.e. due to turbulence the popu- 
lation c(x, t) only saturates locally at a few isolated points. 
The reduction in {Z{t)) tracks in {S{t)), but is much more 
pronounced. 

In Fig. ((T3)) we show Vm{t) computed for the case fi = 
1 and F = 3.. As in Fig. dSJ, we plot Vm * 10 and S{t). 
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Fig. 11. Computation of (S) as a function of F for /i = 1. In 
the insert, we show the prrobability distribution P(S), obtained by 
the numerical simulations, for three different values of F, namely 
F = 0.4,0.8 and F = 2. 




Fig. 12. Computation of the total bacterial mass (Z) (normalized to 
1 at _F = 0) for ^ = 1 (red circles) and of {S)/Smax (green squares) 
as a function of F 

The qualitative behaviour is quite close to what already 
discussed for the case fi = b = 0. The whole picture for ^ = 
1, as obtained by inspection of Fig.s (fTTI) and (fT5|) . supports 
our previous conclusions that, as long as the systems is in 
a quasi-localized phase, the effect of fi in Eq. ^ is almost 
irrelevant. Note that for the system to be in the localized 
phase we must require that both conditions pT|) and (|12p 
must be satisfied. 

According to our interpretation, we expect that for in- 
creasing fi the whole picture does not change provided F is 
increased accordingly. More precisely, we expect that the 
relevant physical parameters are dictated by the ratios in 
Eqs. (fTTjl and p2p . To show that this is indeed the case, we 
show in Fig. (fT4|) the results corresponding to those in Fig. 
(jlip but now with fi ~ 5 instead of /i = 1. 

Two clear features appear in Fig. (|14p . First the qualita- 
tive behavior of (S) with increasing F is similar for /x = 5 



Q 





! ^ s(t) - A.... 


If 





5 10 15 20 25 30 35 40 

t 

Fig. 13. Same as in Fig. ^ for ^ = 1, non zero saturation term and 
F = 3.0. 




Fig. 14. Same as in Fig. (Illl l for /i = 5. 



and ^ = 1. This similarity also applies to the probability 
distribution P{S) shown in the insert of Fig. Second, 
there is a shift of the function (5) p towards large values of 
F, i.e. the localized/extended transition occurs for larger 
values of F with respect to the case fi = 1. This trend is in 
qualitative agreement with the condition ()12|1 . 

To make progress towards a quantitative understanding, 
we would like to use pO)) and ([12]) to predict the shift in the 
localized/extended transition (or crossover) for increasing 
/i. For this purpose, we need a better indicator of this tran- 
sition. So far, we used S* as a measure of localization; large 
values of S mean extended states while small values of S 
imply a more sharply peaked probability distribution. For 
^ = 6 = 0, S" is the "entropy" related to the probability 
distribution P{x, t), solution of eq. ([3]). Thus for /i = 6 = 
we can think of S as the " entropy" and of the diffusion con- 
stant D as the "temperature" of our system. This analogy 
suggests we define a " specific heat" Cs ~ DdS/ dD of our 
system in terms of S and D. After a simple computation 
we get using equations ([8| and ([20]): 
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Fig. 15. (Cs) as a function of F for = 1 (red curve with circles) 
and fi = 5 (green thin curve). The blue line with triangles is (Cs) 
for = 5 plotted against F/b^^"^, for reasons discussed in the text. 

Cs{t)= dxP{x,t)[log{P{x,t))- dxP{x,t)log{P{x,t))Y 
Jq Jo 

(22) 

After allowing a statistically stationary state to develop, 
we then compute the time average (Cs) to characterize 
the "specific heat" of our system for a specific value of 
F. It is now tempting to describe the localized/extended 
changeover associated with ([3]) in terms of the "thermo- 
dynamical" function (Cs). In other words, we would like 
to understand whether a change in the specific heat can 
be used to "measure" the extended/localized transition 
with increasing F. The above analysis can be done also 
for fi > , (when c{x,t) is no longer conserved) by using 
P{x,t) = c{x,t)/Z where the "partition function" Z{t) = 
/p dxc(x,t). 

In Fig. we show {Cg) as a function of F for /i = 1 
(red curve with circles) and /i = 5 (green thin curve). Two 
major features emerge form this figure. First, {Cg) is al- 
most for small F i.e. in the extended case. In the vicinity 
of a critical value F ^ Fc, (Cs) shows a rapid rise to large 
positive values and it stays more or less constant upon in- 
creasing F. The large value of (Cs) reflects enhanced fiuc- 
tuations in log{P{x,t)) (analogous to energy fluctuations 
in equilibrium statistical mechanics) when the population 
is localized. This behavior is in qualitative agreement with 
the notion of phase transition where (within mean field 
theory) the specific heat rises after a transition to an "or- 
dered state". Here, the "ordered state" corresponds to a 
quasilocalized, or sharply peaked probability distribution 
P{x,t). Our numerics cannot, at present, distinguish be- 
tween a rapid crossover and a sharp phase transition. 

The second interesting feature emerging from Fig. (jlSp 
is that Fc, the value of F corresponding to the most rapid 
arise of (Cs) f, depends on /i, as predicted by our theoreti- 
cal considerations. Indeed, as shown just below Eq. we 
expect that Fc ~ (fj,)^^^. To check this prediction, we plot 
in Fig. (|15p a third line (the blue line with triangles) which 
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Fig. 16. Plot of {S) as a function of uq for n = I and F = 2.4. In the 
insert we show the probability distribution P{S) for two particular 
values of uq, namely 1. and 2.2. 

is just {Cg) for /X = 5 plotted against This rescaling 

is aimed at matching the position of the extended/locahzed 
changeover for the same Fc independent of /x. The corre- 
spondence between the two curves in Fig. ()15p confirms our 
prediction. 

Fig. (fT5|) shows that the statistical properties of c{x, t) 
can be interpreted in terms of thermodynamical quantities. 
How far this analogy goes, is left to future research. The 
quantity (Cs) is in any case a sensitive measure of the ex- 
tended/localized transition with increasing F. 

4. Numerical simulations for Uq ^ 0. 

As discussed in Sec. 2, for fixed large F, a mean back- 
ground flow Uq ^ can eventually induce a transition from 
localized to extended configurations of c{x,t). More pre- 
cisely, for large F, i.e. for F large enough to satisfy (fTO]) 
and p^ . the system will spend most of its time in local- 
ized states provided the condition (jlSp is satisfied. Thus, 
for large enough uq we expect a transition from quasi lo- 
calized (i.e. sharply peaked) to extended solutions. In this 
section we study this transition and check the delocaliza- 
tion condition in (fTSl) . 

For this purpose we fix /i = 1 and F ~ 2.4 which, ac- 
cording to our results in the previous section, correspond 
for Uq = to the case where localized states of c dominate. 
As before, we use (S) and of P{S) to characterize the sta- 
tistical properties of c for different values of uq- In Fig. 
we show (S) as a function of uq while in the insert we show 
the probability distribution P{S) for two particular values 
of Uq. For Uq ~ 2 we observe a quite strong increase of (S) 
, a signature of a transition from predominately localized 
to predominately extended states. An interesting feature 
of P{S) for Uq = 2.2 is the long tail towards small values 
of S. This means that, occasionally, the system recovers a 
localized concentration distribution, as if uq = 0. 
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Fig. 17. Plot of the probability distribution P{S) for uq = 1.8, /i = I 
and F = 2.4. In the insert we plot {Cs) as a function of uq, where 
Cb is computed using Eq. I I22II . 

The most striking feature appears near the critical value 
of Uq where the transition a sharp rise in (S')„q occurs. In 
Fig. pT)) we show a two-peaked probability distribution 
P{S) for Uq = 1.8, where the slope of {S)^^ is the largest, 
and in the insert we show {Cs) as a function of uq, where 
Cs is computed using Eq. ((22)) . Let us first discuss the re- 
sult shown in the insert of Fig. p7|) . The specific-heat like 
quantity rises form 0.8 at small uq, shows a bump where 
extended and localized states coexist, and then drops to 0.4 
for uq large. Note that the behavior of (Cg) is different from 
what we observe in Fig. (|15p suggesting a behavior reminis- 
cent of a first order phase transition. We estimate uq ^ 1.8 
as the critical value of uq where the behavior changes more 
rapidly. At uq = 1.8 the probability distribution is clearly 
bimodal, i.e. we can detect the two different phases of the 
system, one characterized by highly localized states and 
the other characterized by extended states. Turbulent fluc- 
tuations drive the system from one state to the other. The 
two maxima in P{S) are suggestive of two different statis- 
tical equilibria of the system. Note that {Cs) is once again 
a good indicator of the transition from predominantly lo- 
calized to predominantly extended states, as discussed in 
the previous section. 

For Uq 7^ 0, a straigthforward generalization of Eq. (|16p 
leads to the following results for the velocity of a maximum 
in c{x, t), 

Vm = fJ,Z / dx{Xm~'X)P{x,t)'^ + u{x , t)P {x , t)dx + UQ 

Jo Jo 

(23) 

One can wonder whether even for uq > 0, the localized 
regime of small S shown in Fig. (fT7|) can be still character- 
ized by Vm ^ 0, thus representing a pinning of the concen- 
tration profile despite the drift velocity uq. This question is 
relevant to understand whether the maxima for small S in 
P{S) shown in Fig. pT)) can be described using ideas devel- 
oped for quasi localized probability distributions in Sees. 
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Fig. 18. Time dependence of S{t) (lower panel) and Vm{t) (up- 
per panel red triangles) for the case of a period mean fiow 
JiO = 1-6 + 0.8cos{2TTt/T) with T = 10, with the same conditions as 
in Fig l|17| l . The line with blue squares in the upper panel represents 
cos(2nt/T) + 1. 




Fig. 19. Contour plot of C{x, t)/Z (the horizontal axis is t while the 
vertical axis is x) for the simulation shown in Fig. I I18I I. 

2 and 3. To answer the above question, we performed a 
numerical simulation with a time-dependent uniform drift 
UQ{t) = 1.6 + 0.8 * cos{2TTt/T) where T = 10.0. Thus 
uq changes periodically in time with an amplitude large 
enough to drive system from one regime to the other. If our 
ideas are reasonable, both S and Vm will become periodic 
functions of time. In particular, as S\ switches from small 
to large values, will go from in the localized regime to 
a large positive value in the extended phase. 

Fig. (jlSp represents a numerical simulation for both S 
and In the upper panel we plot (red line) and the 
periodic function cos{2TTt/T) + 1 (we add an offset of 1 in 
order to make the figure more readable). As one can see, Vm 
indeed flattens out near periodically in time, and increases 
to large positive values in synchrony with the external time- 
dependent drift velocity. In the lower panel, we plot S* as a 
function of time; the graph clearly shows a periodic switch- 
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Fig. 20. The figure shows four snapshots of c{x,t) taken from Fig. 
jigj at times t = 45, 50, 55, 60. At t = 45. and t = 55, when uq = 0.8, 
the population c{x, t) is strongly localized while at t = 50 and t = 60, 
when no = 1.6, c{x,t) is extended. 

ing between the two statistical equilibria. A better under- 
standing of the dynamics can be obtained from Fig. (fT9|) . 
where we show a contour plot of the normalized bacterial 
concentration c{x,t)/Z (the horizontal axis is t while the 
vertical axis is x). Localized states can be observed in the 
vicinity oi t = 45, 55 and t ^ 65 i.e when wq is near its 
smallest value, uq ~ 0.8. Localized states are stationary or 
at most slowly moving whenever uq is small. During the 
period when uq is large, no localization effect can be ob- 
served. Fig. (PO)) we show four snapshots of c(a;, t) taken 
from Fig. HH) at times t 45,50,55,60. At t = 45. and 
t = 55, when uq — 0.8, the population c{x, t) is strongly lo- 
calized while at i = 50 and t — 60, when uq = 1.6, c(a;, t) is 
extended. The reason why fm ^ even for a small Uq > 
is quite simple: according to our analysis in Sec. 2, localized 
states will form near shifted zero velocity points with nega- 
tive slopes even for uq > 0. When uq is large enough, there 
is no point where the whole velocity u{x, t) + uq is close to 
0. Every point in the fluid then moves in a particular direc- 
tion, and the system develops extended states. Fig.s (fTS]) . 
(fTO|) and clearly support this interpretation. 



5. Conclusions 

In this paper we have studied the statistical properties 
of the solution of Eq. ([3|) for a given one dimnesional turbu- 
lent flow u{x, t). Fig. (PT|) illustrates one of the main results 
discussed in this paper: the spatial behavior of the popula- 
tion c{x, t) subjected to a turbulent field. In particular, the 
figure shows four snapshots of c{x, t) taken from a numer- 
ical simulation {D = 0.005, ^ = 1, F — 1.2) at times t = 
60, 65, 70, 75. The population c{x, t) shows strongly peaked 
concentration at time t = 65 and t = 70, while at times 
i = 65 and t = 75, c{x,t) is more extended. The popul- 
tion c{x, t) alternates strongly peaked solutions and more 
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Fig. 21. The figure illustrates one of the main results discussed in 
this paper: the spatial behavior of the population c{x,t) subject to 
a turbulent velocity field. The figure shows four snapshots of c(x, t) 
taken from a numerical simulation (D = 0.005,/^ = I, F = 1.2) at 
times t = 60, 65, 70, 75. The population c{x, t) shows strongly peaked 
concentration at time t = 65 and t = 70, while at times f = 65 
and t = 75, c(x, t) seems to be less peaked. The popultion c{x, t) 
alternates strongly peaked solutions and more extended ones. 

extended ones Our model is sufficiently simple to allow sys- 
tematic investigation without major computational effort. 
From a physical point of view, the model can be interest- 
ing for compressible turbulent flows and whenever the field 
c represents particles (such as the cells of microorganisms) 
whose numbers grow and saturate while diffusing and ad- 
vecting. Our aim in this paper was to understand the sta- 
tistical properties of c{x, t) as a function of the free param- 
eters in the model. We developed in Sec. 2 a simple theoret- 
ical framework. Based on three dimensionless parameters, 
we have identified three conditions which must be satisfied 
for quasi localized solutions of ([3]) to develop, given by Eqs. 

(Uni),!!!]) and nil). 

All numerical simulations have been performed by using 
a grid resolution of iV = 512 points and a Reynolds num- 
ber Re ~ 10^. Increasing the resolution will not change 
the numerical results provided the appropriate rescaling on 
conditions p0)) . p2l) and flS]) are performed, as shown in 
the following argument: let us define Sx the grid spacing, 
i.e. dx = L/N, 77 the Kolmogorov scale and e the mean 
rate of energy dissipation, where 77 = (j^'^/e)^/*. The turbu- 
lent field u{x, t) must be simulated numerically for scales 
smaller than the Kolmogorov scale. In the shell models, this 
implies that the largest value of kn is much larger than I/77. 
The velocity gradient T is of the order of \/Tjv. If the grid 
spacing bx is smaller than 77, no rescaling is needed in us- 
ing the theoretical considerations derived in Sec. 2, namely 
equations pO)) . ([12]) and p8)) . On the other hand, if bx is 
larger than 7/, as in our simulations, the velocity gradient 
goes as 

r ^ + ^ ,l/3(J,)-2/3 (24) 

Thus, by increasing the resolution, i.e. decreasing bx^ we 



Fig. 22. The "specific heat" {Cs) computed for N = 128, D = 0.01 
and /I = 1 (red line with squares), rescaled according to (|24p and 
moi l, and compared with the case fi = 1, D = 0.005 and A'^ = 512 
(green line with circles) discussed in the text. In the insert the same 
quantities are plotted without rescaling. 

increase the velocity gradients and the condition pop may 
not be satisfied unless we change D ov F in an appropriate 
way. As an example of the above argument we show in Fig. 
(|22l) (insert) the value of (Cs) computed for N = 128, D = 
0.01 and /i = 1 (red line with squares) and compared with 
the case, used in the main text, D — 0.005, N ~ 512 and 
/-t = 1 (green line with circles) already discussed in Sec. 3. 
For this particular case, we can superimpose the two curves 
by multiplying F for the iV = 128 case by a factor 0.8 which 
comes from equations ([T0|) and ([24]) . The final result agrees 
quite well with the N — 512 case as shown in the same 
figure. 

Similar considerations apply for a non zero mean fiow ug , 
where the localized/extended transition should occours for 
larger values of uq according to (|18p . Finally, let us mention 
how we can predict the Reynolds number dependence of our 
analysis. According to the Kolmogorov theory, a typical ve- 
locity gradient is F ~ Re^l^ . Therefore, the transition from 
extended to localized solutions predicted by (fTO)) can be ob- 
served provided either D ^ Re~^/'^ or F ^ i?e^/^. Hence, 
by increasing the Reynolds number, the extended/localized 
transition eventually disappears unless the diffusion term 
D is properly rescaled. 

Following the theoretical framework discussed in Sec. 2, 
we introduced a simple way to characterized how well c{x, t) 
is localized in space, namely using the entropy-like function 
(j2ip to illuminate the dynamics of the numerical solutions. 
The time average entropy (S) was used to characterize the 
transition from extended to localized for increasing F and 
from localized to extended solutions for increasing uq. We 
also found it useful to define a " specific heat" Cg by simply 
computing DdS/ dD, where D plays the role of temperature 
in the system. Notice from Eq. (fTU)) that the physics is 
controlled by an effective temperature D/F, i.e. rescaling 
F is equivalent to changing D. 

The analogy between F and some sort of effective tem- 



1 o 




so 100 160 20O 250 300 350 400 450 500 





50 100 160 ZOO 260 300 350 400 450 500 




50 100 150 20O 250 300 350 400 450 500 



Fig. 23. Numerical simulations performed with /i = 1 for 
(7/16)L < a; < (9/16)L and ^ = -1/15. elsewhere for all x. We plot 
the total number of microorganisms Z(t) for three different values 
of F, namely F = 0.25 (upper panel), F = 0.4 (middle panel) and 
F = 0.5 (lower panel). 

perature suggests that the rapid rise in the time average 
(Cs), observed in Fig. (fT5|) near a characteristic value Fc, 
might indicate a critical "temperature" or diffusion con- 
stant Dc- Fig. (fT5|) highlights the rapid changes in {Cs)f 
from extended to localized states in the system. It will be 
interesting to study the behaviour observed in Fig. (|I5|1 
from a thermodynamic point of view. As predicted by pre- 
vious analytical studies [6] [7], with time- independent ve- 
locity field, a transition from localized to extended states 
has been observed by increasing uq. The interesting fea- 
ture is that near this transition, the system shows a clear 
bimodality in its dynamics, at least in the probability dis- 
tribution P{S), more indicative of a first order transition. 

We are not able at this stage, to predict the shape of the 
probability distribution of P{S) as a function of external 
parameters such as D,F,fi and uq. It would be valuable to 
understand better when a quenched approximation (time- 
independent accumulation point in u{x,t)) is reasonably 
good for our system, especially in regimes where microor- 
ganism populations are nearly localized. The reason why 
a quenched approximation may work is that the localized 
regime is quasi-static, in the sense that the solution c(x, t) 
follows the slow dynamics of accumulation points where 
u{x, t) ~ Q with a large negative slope. A complete discus- 
sion of the validity of quenched approximation and analytic 
computations of P{S) is a matter for future research. 

So far we have discussed the case of ^ constant and pos- 
itive. In some applications (both in biology and in physics) 
one may be interested to discuss /i with some non triv- 
ial space dependence. An interesting case, generalizing the 
work in [6] , [7] and [12] is provided by the equation. 



dtc + d^{Uc) = Ddlc + ^i{x)c - bc^, 



(25) 



with a turbulent convecting velocity field U{x,t) = uq + 
u{x, t) and where /i is positive on a small fraction of the 
whole domain and negative elsewhere. In this case, referred 



to as the "oasis" , one would like to determine when (c(a-, t)) 
can be significantly different from zero, i.e. when do the 
populations on an island or oasis survive when buffeted by 
the turbulent flows engendered by , say, a major storm (see 
[8^ for a treatment of space-independent random convec- 
tion). A qualitative prediction for (c(x, t)) results from the 
following argument: the extended and/or localization be- 
haviour of c depends on the ratio defined in Eq. ((12]). For 
small F (i.e. small F) the solution must be extended and 
therefore one can predict that c is significantly different 
from zero everywhere, wherever /.t > 0. On the other hand, 
for large F. c becomes localized. The probability for c to be 
localized in one point or another is uniform on the whole 
domain. Thus, if the region where > is significantly 
smaller then the region where /.t < 0, c should approach to 
zero for long enough time. 

In Fig. (|23p . we show a numerical simulation performed 
for a oasis centered on x = i/2, performed with /i = 1 
for (7/16)L < X < (9/16)i and = -1/15. elsewhere 
for all X. Thus the spatial average of the growth rate is 
dxfi{x) = 1/15.. As a measure of c, we plot its spa- 
tial integral Z{t) as a function of time. The numerical sim- 
ulations have been done with D = 0.005, N = 512 and 
Uq = 0. In Fig. ([23]) we show three different values of F. Let 
us recall that, when fJ- = 1 everywhere, as is now the case 
for the oasis, for F > 0.25, the system exhibits a transition 
from extended states to localized states, as illustrated in 
Fig. ([TS]) . As one can see, for F > 0.4 the population tends 
to crash as predicted by our simple arguments. It is how- 
ever interesting to observe that the dynamics of c is not at 
all trivial. For F = 0.4 and F ~ 0.5, c seems to almost die 
and then recovers. Of course, our continuum equations ne- 
glect the discreteness of the population. At very low pop- 
ulations densities, a reference volume can contain a frac- 
tional number of organisms and cxtintion events are artifi- 
cially supressed. Fig. ([^ neverthclesst suggests an inter- 
esting feature of Eq. ([25]) with space dependent /i, worth 
investigating in the future. 

Another interesting question tjhat deserves more de- 
tailed studies is the case two competing species with densi- 
ties ci{x, t) and C2{x, t). To illustrate the problem, consider 
the coupled equations: 

dtci + d^{Uci) = Ddlci + A*ici(l - ci) - /i2CiC2 (26) 
dtC2 + d^{Uc2) = Ddlc2 + A*2C2(1 - C2) - /iiCiC2 (27) 

where 112 > Mi and < (5/i = /i2 ^ /^i <C /xi. In this 
simplified model, Eqs. (1261) and ([27]) describe the dynamics 
of two populations in which a "mutant" density C2 can 
out compete a wild type density ci. In particular, upon 
specializing to one dimension and denoting Csum{x,t) = 
ci(x, t) + C2{x, t), from Eq.s (|26I27|) we obtain: 

dtCsurn + dx{U Csum) = Ddlcsum + (1 - Csu7n){p-lCl + H2C2) 

(28) 

Eq. (|28p shows that, for U ~ 0, Cs,„„ = 1 is an invariant 
subset, i.e., if at i = 0, Cgum = 1, then Csum = 1 for any 
t. The system has two stationary solutions, namely ci = 
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Fig. 24. Numerical simulation of equations 1 1261 1. 11271 with = 1. 
and fi2 = 1.05. We show the quantities Z1/Z2 (open circles) for 
the case with no turbulence. Note that Z1/Z2 decays to zero as 
ea:p(— (/J2— (green line) . When turbulence is acting, the dynam- 
ics becomes more intermitent as shown by the behavior of (Zi)/ {Z2) 
(close red circles) and and {Z1/Z2) (solid triangles). The symbol (...) 
means averaging over ensemble. In the insert, we show the proba- 
biltiy distribution of log[Zi(t) /Z2{t)], at t = 20, which is well fitted 
by a guassian behavior. 

1,C2 ~ which is unstable, and ci = 0,C2 ~ 1 which is 
stable. For U = 0, any initial conditions is attracted to the 
stable solution. It is easy to check that the asymptotic time 
dependences in this subspace are of Zi ^ exp{—6^t) and 
Z2 ^ I — exp{—S^t), where Zi = J dxci{x, t). 

In Fig. (|24p we show the result of two different numeri- 
cal simulations of Eqs. (|26I27|) with /ii = 1 and yU2 = 1.05. 
The solutions have been obtained by using periodic bound- 
ary conditions, L = I, D — 0.005 and a numerical resolu- 
tion of 512 grid points. The open circles represent the be- 
havior of log{Zi{t)/Z2{t)) for [/ = 0. As predicted by our 
simple analysis, Z1/Z2 decays quite rapidly towards as 
exp{—5iJLt) (dashed green line in Fig. ([24|) . Note that — 
0.05, corresponding to a characterstic time 1/Sfi ^ 20. 

For U ^ 0, however, the time behavior is quite different. 
In particular, we choose uq = and allow convection by a 
strong turbulent field with F = 0.8. In Fig. ([21]), the red 
circles refer to log{{Zi{t)) / {Z2{t))) while the blue triangles 
refer to log{{Zi{t) / {Z2{t)) . The symbol (...) is the ensemble 
average over 100 realizations of the turbulent field, with the 
same initial conditions 

ci(x,t = 0) = 1 C2(x,t = 0) = /or Q<x<^ (29) 
Ci(.T,< = 0) = C2(a;,< = 0) = 1 /or ^<x<L (30) 

While the asymptotic states are still the same as for the 
case F = (the stability of the stationary solutions does 
not change) , the population ci decays on a time scale longer 
than the F = one (i.e. \/5ii). The rather large differ- 
ence between {Zi{t)) / {Z2{t)) and {Zi{t) / Z2{t)) is due to 
strong fluctuations in the ensemble. To highlight these fluc- 
tuations, we show in the insert of Fig. ([M)) the proba- 



bility distribution P{R) of the logarithmic ratio R{t) = 
log{Zi{t) / Z2{t)) computed at t = 20, which is well fitted 
by a gaussian distribution with a rather large variance. 
This implies that the ratio Z1/Z2 is a strongly intermit- 
tent quantity. To explain such a strong intremittency, note 
that the initial time behavior of the system strongly de- 
pends whether one of the two populations is spatially ex- 
tended while the other sharply peaked. When the popu- 
lation ci(a;,t) is extended while C2{x,t) is sharply peaked, 
the ratio Zi/ Z2 becomes initially quite large. On the other 
hand, when C2{x,t) is extended and ci(x, t) sharply peaked, 
Z1/Z2 is very small. For long enough times, the two popu- 
lations become correlated in space (by clustering and com- 
peting at the same accumulation points of u{x, t)) and the 
ratio Z1/Z2 eventually decays according to the expected 
behavior exp{^S^,t). Note that the characteristic turbu- 
lent mixing times in our simulations are much longer that 
the characteristic doubling times of the microorganisms, 
~ 1/fJ-i and ~ 1/M2- This is the opposite of the situa- 
tion in many microbiology laboratories, where organisms 
in test tubes arc routinely mixed at a rapid rate overnight 
at Reynolds numbers of the order lO'^. The situation stud- 
ied here can, however, arise for microorganisms subject to 
turbulence in the ocean. 

We close with comments on generalization to more that 
one dimension. When "turbulent velocity field" u{x, .t) can 
be represented as V5'(2:, t), with a suitable 4", most of the 
results discussed in this paper should be valid. However, in 
a real turbulent flow in higher dimensions, whether com- 
pressible or incompressible, the velocity field is not irrota- 
tional. For a real turbulent flow, we believe the localization 
discussed here will be reflected in a reduction of the space 
dimensions in the support of c{x^t). For instance, in two 
dimension, we expect that c(a;,t) will become large on a 
onc-dimensional filament while in three dimension c{x,t) 
localizes on a two dimensional surface. For a review of re- 
lated effects for biological organisms in oceanic flows at 
moderate Reynolds numbers see |17] 

Following the multifractal language, there may be a full 
spectrum of dimensions which may characterize the statis- 
tical properties of localized states. It remains to be seen 
whether a sharp crossover (or an actual phase transition) 
similar to what has been shown in Sees. 3 and 4, will be 
observed in more than one dimension. 
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